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^\ , The effect of a bulk marginal operator on boundary critical phenomena in two space-time di- 

■ mensions is considered. The particular case of an open S=l/2 antiferromagnetic Heisenberg chain, 
corresponding to a Wess-Zumino-Witten non-linear a model, is solved. In this case, the needed 
renormalization group coefficient is associated with a novel operator product expansion in which 
three operators approach the same point. Resulting logarithmic corrections occurring in finite size 
calculations and nuclear magnetic resonance experiments are discussed. 

CO ■ PACS number 75.10.Jm 

D 1 I. INTRODUCTION 

h. 

Marginally irrelevant operators in two-dimensional conformal fieldm theory lead to logarithmic corrections to scaling 
behaviour. Because the corresponding coupling constant, g(l), renormalizes to very slowly as 1/lni where I is a 
characteristic length or energy scale, logarithmic corrections occur to virtually all quantities which can be measured 
experimentally or simulated numerically. These create grave difficulties in obtaining agreement between analytical 
' theory and numerical simulations or experiment. The particular case of the S=l/2 Heisenberg antiferromagnqtig. 
chain has been discussed extensively. The correlation function, initially predicted to decay as 1/r, instead decays asBcl 
(In r) 1 / 2 /r. To make matters worse, the corrections to this result are only suppressed by additional powers of 1/lnr 
and are highly sensitive to finite-size effects. Similarly, the energy gap between the ground state and the first excited 
(triplet) state behaves asu: 
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■ where v is the spin-wave velocity. At very long lengths, 
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giving an additive logarithmic correction to the finite-size energy gap. It is generally very hard to actually observe this 
logarithmic behaviour unless chains of length several thousand can be studied. Fortunately this is possible for Bethe 
ansatz integrable models like t he S =l/2 Heisenberg chain. For shorter chains it is generally better to regard g(l) as a 



free parameter. Because Eq. (1.1) has been generalized to many other energy levels, all of which receive corrections 
linear in g, this still has considerable predictive power. Indeed, fitting to expressions of this sort, which effectively 
subtracts off the leading logarithmic.correction to scaling, provides a practical method for numerically determining the 
J-J ' universality class of a HamiltonianB There are analogous finite temperature corrections (for infinite length systems). 
, These include an additive 1/ InT correction to the susceptibilitya and a multiplicative (InT) 1 / 2 correction to the spin 
relaxation rate, 1/Ti.Q Although experimental data has been fit to these forms with apparent successB'Lj this is a very 
difficult exercise due to the slow variation and the presence of various other types of corrections in any real material. 



Another subject of current experimental and theoretical interest is the general area of quantum impurity problems 
(QIP's). In the context of S=l/2 Heisenberg antiferromagnetic chains, such a problem occurs for a semi-infinite chain, 
associated with the dynamics at the chain end.EI This model can be realized experimentally by dilute substitution 
of the magnetic ion by a non-magnetic one, eg. Zn substitution for Cu. Related QIP's involve the Kondo problem 
and tunneling through a single impurity in a quantum wire (or quantum Hall effect edge states). The general 
renormalization group (RG) treatment of these problems gives fixed paints corresponding to conformally invariant 
boundary conditions imposed on a given bulk conformal field theory.B In general, in such a theory, the effective 
Hamiltonian contains both bulk and boundary operators. The bulk terms contain integrals over the half-line whereas 
the boundary terms occur at the impurity location, x — 0. Bulk behaviour is unaffected by boundary dynamics, 
although it must be appreciated that the decay of Green's functions away from the boundary is itself part of the 
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boundary critical phenomena. Time correlation functions at the boundary also involve exponents which characterize 
the boundary condition, as does the finite size spectrum with non-periodic boundary conditions. Boundary interactions 
cannot affect the renormalization of bulk coupling constants. In most treatments of these problems so far, any 
renormalization of boundary interactions by bulk interactions has also been ignored. The justification for this is that 
the bulk system has been assumed to be at a bulk RG fixed point. Any bulk operators present, apart from the fixed 
point Hamiltonian itself, are irrelevant and can be ignored at low energies. Thus, at least in principle, crossover 
between boundary fixed points can be treated independently of bulk renormalization. Strictly speaking, this is only 
justified when the energy scales associated with the boundary renormalization are much smaller than those associated 
with the bulk irrelevant couplings. This approximation is particularly bad when there are marginally irrelevant bulk 
interactions present since they renormalize to logarithmically slowly. When marginally irrelevant bulk operators 
are present we should expect logarithmic corrections to exponents and finite-size scaling. However, the detailed form 
of these corrections is characteristic of the boundary condition and is not simply related to the log corrections in the 
bulk theory, nor to the finite-size scaling with periodic boundary conditions. 

It is the purpose of this note to consider these logarithmic corrections to boundary critical exponents and finite-size 
scaling with non-periodic boundary conditions arising from a marginal bulk operator. In the next section we calculate 
general formulas for the logarithmic corrections using conformal field theory. In Section 3 we compare these formulas 
to results on the finite size spectrum for an S=l/2 chain with open boundary conditions, for lengths up to 2,000, 
obtained from the Bethe ansatz. In the final section we comment on the log corrections to correlation functions and 
the NMR relaxation rate. 



II. CONFORMAL FIELD THEORY RESULTS 

As observed by CardyS such logarithmic corrections generally result from corrections to the anomalous dimensions 
of the various operators, (f> n , which are linear in g. The associated coupling constants, «„, obey an renormalization 
group (RG) equation: 



du n /d\nL= [2- 7n ]<7, (2.1) 



where: 



7n = x n + 2nb n g + ... (2.2) 

Here the . . . represents terms of higher order in g or in other irrelevant operators. Taking into account the fact that g 
itself also renormalizes leads to predictions of the various logarithmic corrections. In particular, the finite-size states 
are in one to one correspondence with the operators and their energy gaps are proportional to j n . The coefficients, 
b n can be conveniently determined from the operator product expansion (OPE) of the operator <f> n with the marginal 
operator, <f>: 

4>{z)K{z') ^ ^j^- (2.3) 
A calculation of any Green's function involving <j) n , to the first order in <?, encounters a logarithmic ultraviolet 



divergence upon integrating z near z' . This implies the correction to the anomalous dimension in Eq. (2.2). The 
finite size spectrum is given by: 

E n -E ~— [x n + 2Trb n g(l)} . (2.4) 

The scaling dimension x n of the corresponding operator is simply corrected by the anomalous dimension term of first 
order in g which is replaced by the effective coupling at scale I. 

In the presence of a boundary, this calculation takes a rather unfamiliar turn because the local marginal bulk 
operator in general becomes bilocal in the presence of a boundary condition. This follows from Cardy's general 
approach to conformally invariant boundary conditions, which are always assumed to obey: 

T L (t,0)=T R (t,0), (2.5) 

where r are the left and right moving terms in the energy momentum tensor. Since Tl r is a function of (t — x) 
[{t + x)] only, it then follows that we may regard the right-movers on the original physical space, x > as the 
continuation of the left-movers to the negative axis, 
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T R (x)=T L (-x), (x>0). 



(2.6) 



This observation allows the Hamiltonian to be written in terms of left movers only, but defined on the entire real 
line. In particular, it implies that a generic local bulk operator, which can be factorized into its left-moving and 
right-moving parts, 0\(t — x) and 2 R (t + x) respectively, becomes bilocal: 

0{t, x) = 0\{t, x)0 2 R {t, x) -» Ol(t, x)0 2 L (t, -x). (2.7) 

In particular, the bulk marginal operator becomes bilocal, introducing a novel complication in calculating its effects 
perturbatively. Boundary operators are also drawn from the left-moving sector only. We may determine the correction 
to the anomalous dimension of an arbitrary boundary operator, <E>, to first order in g, from the 3-point function 
< $0$^ >, where O is the marginal bulk operator. However, this 3-point function mu st b e calculated in the presence 
of the boundary condition, upon which it depends. Furthermore, we see from Eq. ( |2.7| ) that this 3-point function 
effectively becomes a 4-point function of left-moving operators. 

In general, some data about the boundary condition will be needed to calculate this 4-point function. As will be 
seen from the example considered below, it is sufficient to know the OPE of the chiral part of the bulk marginal 
operator, O(z) with general boundary operators. We consider a special case here, of some importance, for which this 
OPE can be readily calculated. 

Let us consider the problem of an S=l/2 Heisenberg chain with an open boundary condition. The bosonized form of 
this model, applicable at low energies, is the k — 1 Wess-Zumino-Witten non-linear a model. The marginal operator 
is quadratic in the chiral spin densities, Jl,r, and is written: 

H = H - 9 (8t: 2 /V3)J l ■ Jr. (2.8) 
Here the spin densities are normalized as: 

< J£(z)j|(0) >= (2.9) 

The factor of 8n 2 /v3 is inserted so that the operator multiplied by g in the Hamiltonian has a unit-normalized 2-point 
function, following the convention of Cardyiil The needed OPE's in this case follow from the basic one of the WZW 
model: 

Jl(z)^z') = -^1- + ... (2.10) 

2n{z — z') 

Here the general Virasoro primary operator, 4>, may transform under an arbitrary representation of SU(2) l x SU (2)r. 
The finite dimensional matrices, are simply the representation of SU(2) under which ip transforms. Explicitly, if 
<p transforms under the irreducible spin S representation then there will be a multiplet of (2S + 1) operators, 4>a and 
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Thus, 



J L {z) ■ J R (z)^(z') = -^-^d>{z) + ... (2.12) 



and we see that the coefficient b n define in Eq. (2.3) takes the value: 



K = - 2 -^. (2.13) 
Now consider a semi-infinite chain with a free boundary condition at the origin. This bound ary condition was 



treated using bosonization in Ref. [ |8J. The boundary condition not only is consistent with Eq. (2.6) but also with 
its Kac-Moody generalization: 

Jr{x) = J L (-x). (2.14) 



Eq. (2.6) and (2.14) imply that the correlation function of T(t, x) and J(i, x) are simply the chiral correlation functions 



in the free WZW model. The boundary has no affect on them apart from identifying left with right. The operator 
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content and finite size spectrum are drawn from conformal towers of the left-moving Kac- Moody algebra only. It was 
shown in [ |§| that, for an even length chain, the spectrum consists of the identity conformal tower only. In particular, 
the spin operator at the boundary becomes the (left-moving) spin density, Jl (t, 0) of the WZW model with correlation 
function, oc 1/t 2 . The finite-size spectrum, with free boundary conditions at both ends of a chain of length I is given 
by: 

E-E = —x L , (2.15) 
where xl is the scaling dimension of t he corresponding (left-moving) operator. This differs from the formula with 



periodic boundary conditions, Eq. (2.4), by the replacement of I by 21 corresponding to doubling the system size due 
to the identification of Eq. (2.6) at x = and x = I and to the replacement of the scaling dimension x = xl + xr by 
a left-moving scaling dimension xl only. 

Let us now consider the logarithmic corrections to this formula for open boundary conditions. The problem again 
reduces to finding the correction to the anomalous dimension of a given operator, <p n , to first order in g where <f> n is 
now a boundary operator. We might again attempt to obtain this from an OPE but we must deal with the fact that 



the marginal operator is now bilocal. The obvious generalization of Eq. (2.12) is 



J L (z) ■ J L (z*)m « ^^0- (2.16) 

However, this is not a conventional OPE because we are bringing three operators to the same point rather than just 
two. However, as will be argued below, it is nonetheless correct. It is now easy to calculate the correction to the 
anomalous dimension. A logarithmic ultraviolet divergence is again encountered for any Green's function involving 
the operator <fr, to first order in g from integrating over z near 0. However this is reduced by a factor of two relative 
to the bulk case because z is only integrated over the half-plane rather than the entire plane. A convenient ultraviolet 
cut-off is given by restricting the z integral to \z\ > a. In the bulk case the excluded region is a circle of radius 
a around the origin but in the boundary case it is only a se mi-c ircle of radius a. Our conclusion is thus that the 



anomalous dimension of a boundary operator is given by Eq. (2.2) but with the coefficient b n now given by 



b n = -^L. (2.17) 



This differs from the bulk formula Eq. (2.13) only by the identification of Sr with Sl and by the extra factor of 1/2 
arising from the different integration region. 

To complete our derivation we just need to justify the rather unorthodox 3-operator OPE occurring in Eq. (2.16). 
The validity of this formula can be understood by considering the more general connected 4-point function: 

G AB =< J L ( Zl ) ■ J L (z 2 )0»0 B (r) > connected • (2.18) 

We normalize the Virasoro primary boundary operator, <j) A so that its two point function is given by: 

SjAB 

< ^(°) At) >= (2.19) 
We wish to consider the short-distance singularity in G AB when Z\,Z2 —* 0. We claim that this is given by: 

G AB -|^7^- (2-20) 

47T z Zi2:2 ( — T) 

The correctness of this result can be seen by considering the three different limits \z\\ << \zq\, \z\ — Za| and \z%\ << 
\zi\, \zi — #2 1 and \z% — z%\ << \zi\,\z2\. In the first case, we can obtain the leading singularity by using the OPE 
of Jl{z\) with (0) and then the OPE of the result with Jl{zi)- This gives Eq. (2.2C). The same singularity is 



obtained in the second case. In the third case there should be no singularity of the form 1 / [z\ — z<z) because there is no 
singular term in the OPE Jl{zi) ■ Jl{z2) apart from the trivial one which doesn't contribute to the connected Green's 
function. These considerations uniquely fix all singularities in G AB at zi,Z2 — > 0. Note that the crucial property of 
the boundary condition that is being used is that the OPE of the spin density operators </£ with arbitr ary ( Virasoro 



primary) boundary operators has the same form as in the bulk. Now letting z\ = z, Z2 = z* gives Eq. (2.16) 
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In the particular case where the <pA are the spin-density operators, J£, we have calculated G ab exactly and verified 
the form of the singularity. In this case we find: 



< Jl(zi)- J L (z2)J a L (0)J b L (r) > = 



gab 



(27r) 4 ziz 2 (T-z 1 )(r-z 2 )' 



(2.21) 



We see from Eq. (2.9) that the unit normalized operator is (j) a = 2-K^/2J a . Also using the fact that J has 5j, = 1 and 
therefore Sl • Sl — Sl{Sl + 1) = 2, we see that, in the limit Zi — > 0, Eq. (|2.21|) agrees with Eq. (|2.20|) . 



Now let us consider the finite size spectrum, examining the lowest energy excited state of spin S. This is given by 
Eq. (2.15) except tha t th e dimension of the (left-movi ng) fi eld, xl must be replaced by the anomalous dimension, j n . 
This is given by Eq. (2.2) with b n now given by Eq. (2.17). Thus we obtain: 



E 



open 



„open 

tin ~ — T 



7TV 

T 



s 2 



2nS(S+l)g(l) 



V3 



For exponentially long chains we may use the asymptotic form of g(l): g(l) — > "y/3/ (47r In i) , giving: 
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(2.22) 



(2.23) 



It is interesting to compare Eq. ( |2.23 ) to the corresponding result for periodic boundary conditions.! If we again 
consider the lowest energy state of given spin S, this has Sl = Sr = 5/2 and hence: 
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(2.24) 



The l/l terms are the same for open and periodic boundary conditions but the 1/lnl terms are not. (We note that 
the logarithmic corrections for open boundary conditions were assumed to be same as the ones for periodic boundary 
conditions in [ |icj )1 . 

We have also calculated the logarithmic correction to the ground state energy for open boundary conditions. 
Ignoring the irrelevant operator, the ground state energy for any one dimensional Hamiltonian which renormalizes to 
a co nform al field theory, defined on an interval of length / with generic boundary conditions at and I consistent with 
Eq. (U) is: 



E (l) = en,l + ei - (tt«/24/)c, (2.25) 

where c is the central charge.0Note that the coefficient of l/l is 1/4 times the value for periodic boundary conditions. 
Also note that an additional non- universal surface energy, e± appears when the boundary conditions are non-periodic. 
Logarithmic corrections to this formula can be calculated by doing perturbation theory in the marginally irrelevant 
coupling constant, g(l) just as in the periodic case. One finds that the correction of 0(1/1) is universal. This must 
be separated from various non- universal corrections to eo and e\. Once the correction to the l/l term of leading 
order in g is calculated, g may be replaced by the effective coupling constant at scale I, g(l), resulting in logarithmic 
corrections. In the case of periodic boundary conditions, this leading correction was found to be 0(g 3 ). By contrast, 
in the case of open boundary conditions we find that it is Q(g 2 ) 



Let us first consider the correction of O(g). From Eq. (2.8) this gives a correction to the ground state energy: 



8E, 



(i) 



V3 Jo 



dx < Jl(x) ■ Jl(—x) > . 



(2.26) 



This Green's function is given in Eq. (2.9) for the case I — > oo. We may obtain the Green's function for finite length 
by a conformal transformation: 



(r' + ix') = e (^/0(T+«) ) 



(2.27) 



giving: 



< Jl{x)J b L {-x) >-- 



87T 2 



(2.28) 
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The integral of Eq. ( [2.26 ) is ultraviolet divergent both at x = and x — I. We may insert an ultraviolet cut off on 
the integration region, x>a, I — x>a where a is of order the lattice spacing. This gives the ground state energy 
correction of first order in g: 



(2.29) 



Now Taylor expanding in powers of a/ 1, we see that we obtain a cut off dependent contribution to the non- universal 
surface energy, e\ in Eq. ( 2.25| ) together with corrections of 0(l/l 2 ): 



SE, 



(i) 



-29^3 



0{a/l 2 ). 



(2.30) 



Importantly, there is no term of 0(1/1). 

We now push this calculation to second order in g. This term is given by: 



6EP = 



V3 



dr dxx dx 2 < Jl{t,xi) • Jl(t, -XijJ L {Q, x 2 ) • Jl{0, -x 2 ) > 



Let us first evaluate this expression in the limit I — > 00. Using Eq. ( 2.21 ), we obtain: 

2 



SE, 



(2) 



8tt 2 

7T 



(27T) 



d.T 



dx\dx 2 



o [t 2 + ( Xl ~x 2 ) 2 ][t 2 + (jbi +x 2 ) 2 



(2.31) 



(2.32) 



Note that we have inserted a factor of 2 here because equal contributions arise from x near and x near /. The 
x-integrals can be done exactly and are ultraviolet finite, for non-zero r: 



SE, 



(2) 



g 2 TT 2 r dr 



(2.33) 



The r integral is ultraviolet divergent. We cut off the integral, |t| > To, where To is of 0(a/v). This gives the second 
order ground state energy correction: 



SE, 



(2) 



TO 



(2.34) 



another cut-off dependent contribution to the surface energy, e% in Eq. ( |2.25 ). To obtain SE^ at finite I, we again 
use the conformal transformation of Eq. ( 2.27 ) to obtain: 



SE^ 



dr 



dx\dx 2 {K /2l) 4 



I sin(7r/20[(xi - x 2 ) 2 + ir}\ 2 \ sin(7r/2Z) \{xi + x 2 ) 2 + ir}\< 



(2.35) 



The Xi integrals are again finite for non-zero t. We again cut off the r integral at |t| > To. Noting that the integrand 
is symmetric under X\ — > —X\ or x 2 — > —x 2 and also t — > — t, it is convenient to extend the Xi integrals from —I to / 
and reduce the t integral from To to 00. This introduces a net factor of 1/2. It is now convenient to change variables 
to: 



and 



u = ttt/1. 



(2.36) 
(2.37) 



The complex variables, Zj are integrated around the unit circle. This expression now becomes: 



SE^ = —^-j— I du I dz\ I dz 2 



C Jc z 2 (zi - z 2 e u ){zi - z 2 e u )(zi - z 2 e u )(z 1 - z 2 e ") ' 



(2.38) 



where uo = itt /1 and C denotes the unit circle integration contour. The z\ integral can now be done by the standard 

4' 



contour integration method, with contributions from the poles at z^ x e The result is: 
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5E^ = -Wj 



du coth 
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Also doing the Z2 integral by contour methods gives: 

5E& = -*V(*/0 
Finally, performing this elementary integration gives: 

6E® = -2ttVW0- 



du 



u sinh (u) 



,2uo _ \ 



T V , 2 2^ 

+7T y 

T I 



O(t /1 2 



(2.39) 



(2.40) 



(2.41) 



We have recovered the same cut off dependent to the surface energy, e\ as in Eq. ( 2.34 ). More importantly, we have 
also obtained a term of 0(1/1) which is cut off independent and therefore is expected to be universal. Thus, we obtain 
the log correction to the ground state energy with open boundary conditions: 
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For exponentially large I we may use Eq. (1.2) to write 
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(2.42) 



(2.43) 



As usual, the corrections to these formulas are only down by additional powers of g(l) , that is 1/ In I. The corresponding 
formula for periodic boundary conditions is: 
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(2.44) 



Among other differences, note that the log corrections decrease the apparent value of c for open boundary conditions 
but increase it for periodic boundary conditions. The fact that the correction to c goes like l/(ln/) 2 was obtained 
from the Bcthe ansatz in Ref. [ [l^| , although the coefficient was not obtained. 



III. NUMERICAL RESULTS ON FINITE SIZE SPECTRUM 



One application of the above results for boundary critical phenomena is to the numerical study of finite size 
scaling. We extract estimates of g(l) defined in Eq.( [2.22| ) and Eq.( |2.42 ) from the energies of finite size spin 1/2 
antiferromagnetic Heisenberg open chains. The Hamiltonian is 



i+l) 



(3.1) 



where the Si's are spin 1/2 operators. The Bethe ansatz equationstJ for the Hamiltonian of Eq. (3.1) are 



A fc + i/2 
A fc - i/2 



n 



A k - A j 



+ i A fe + A., + i 



A k -A«-i A fe + A, 



(3.2) 



where I is the number of sites in the open chain. The roots can be numerically calculated@ The number of roots, M 
determines the total S z co mponent through the relation S z = L/2 — M . In terms of the solutions, Afc, to the Bcthe 
ansatz equations, Eq. (3.2), the energy is given by 



E = 



l-l 



1 M 1 

1 V — 

2 ^ A\ + 1 

fc=i k 



/4' 



(3.3) 



The surface energy ei = (tt — 1 — 2 ln2)/4 for the two ends of open chain can be exactly obtained. The rapidities 
{Afc,fc = 1,1/2} solving the Bethe ansatz equations Eq. ([T5) for the ground state of the open chain are all bigger 
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than zero. Let's order A/, so that: Ai > A2 > . . . > A;/ 2 -i > A//9 > 0. We can construct another set of rapidities 

{A' k ,k = i,iy. 

A^-A,-, A£ +1 _,. = A i , for j = l to 1/2, (3.4) 
which solves exactly the following Bethe ansatz equations for a periodic chain of length 21 + 1: 

'k' k + i/2\ 2l+1 _ A' fe - a;. 



M-i/2) 11 K -A', 



The energy for the periodic chain is given by 

, 21+1 1^ 1 , s 

E=— ; -V^ . 3.6 



fe=i 



In the logarithmic form of the Bethe ansatz equations for the periodic chain, the set {A^, k = 1,1} corresponds to a 
set of integers {I' k , k — 1,1}: 

I] = j- (1/2+1), = 3i for j = l to 1/2, 

i.e., -1/2, -1/2 + 1, ...-2, -1, 1, 2, ...1/2-1, 1/2. (3.7) 

The groundstate of the open chain corresponds to an excited state of the periodic chain jsnth a hole exactly at I L = 
in the connected integer set {/,, i = 1, 1 + 1} for the ground state of the periodic chain. O While the groundstates of 
the periodic chain of odd length have spin S z — ±1/2 and momentum ±7r/2, the state of the odd length periodic 
chain corresponding to the groundstate of the even length open chain has S z — —1/2 and momentum 0. From its 
momentum, and from the fact that it has one hole we expect its energy, in the large I limit, to be the groundstate 
energy of the periodic chain plus the excitation energy for a magnon of momentum tt/2: E' = eo(2l + 1) + v s sin(fc), 
v s = tt/2, and eo = 1/4 — In 2 for the spin 1/2 chain. Comparing the expressions for E and E' , we can eliminate the 
summation over A and obtain E = e$l + e\ for the ground state of the open chain with e\ = (w — 1 — 21n2)/4. We 
have checked this result from our numerical solution of the Bethe ansatz equations for finite I obtaining agreement 
to at leastiLdecimal places. This result was derived earlier from the Bethe ansatz equations by a somewhat different 
methodMliile 

To test the CFT predictions, we extract three estimates of g(l) using Eq.( 2.22j ) for the S = 1 a nd S — 2 excited 



states and Eq.( 2.42j ) for the groundstate for chains up to two thousands sites by solving Eq. ( |3.2[ ) numerically. We 



draw the three g(l) completely determined by the energies of these three states, respectively, in Fig. [|. The reason 
that these three estimates of g(l) don't agree exactly is because of the various corrections to these energies of higher 
order in g(l). However, at large I these estimates should converge since g(l) — > 0. We see that the coupling constants 
indeed collapse into one value which approaches 0, as the chain length increases thus verifying the CFT predictions. 
Previous numerical studies for periodic boundary conditions have verified finite size scaling obtained by conformal 
field theory&BEl We then compare the g(l) obtained from the open boundary conditions with the one from periodic 
boundary conditions in Fig. |^. We redraw the ground state g{l) for open boundary conditions and the average g(l) 
given by the ground state, the singlet excitation, and the triplet excitation for periodic boundary conditions. These 
data for periodic boundary conditions were obtained in Ref. [ [l6| . The two g(l) approach each other in the large length 
limit. The one-loop renormalization group prediction for g(l) given in Ref. [ |^], g(l) = go{lo)/[l + itbg^lo) ln(//Zo)], 
with b = 4/\/3 is also drawn in Fig. |^. We use the average of the g(l) for periodic and open chains at Iq — 2048 to 
fix go. The one- loop renormalization group prediction fits g(l) at large I. So we see that the logarithmic corrections 
for Heisenberg chains have been successfully predicted by conformal field theory. 
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FIG. 1. g(l) calculated by Eq.(2.22) and Eq.( p.42| ) from Bethe ansatz energies of spin 1/2 antiferromagnetic Heisenberg 
open chain (OBC). The three g(l) are obtained from energies of ground state, total spin S = 1 excited state and 5* = 2 excited 
state, respectively. 




FIG. 2. The g(l) calculated from ground state of open chain (OBC) and the average g (I) for periodic^ chain (PBG1 
for spin 1/2 antiferromagnetic Heisenberg model. The full line is the one-loop renormalization group (RG) prediction^, 
g(l) — go(lo)/[l + nbgo(lo)ln(l/lo)] with b — 4/y3. go{lo) is determined by the average of open chains ground state's g(l) 
and the periodic chain's average g(l) at the chain length lo — 2048. 



IV. CORRELATION FUNCTIONS AND 1/T] 



Another application of the anomalous dimension of boundary operators is to Green's functions for a semi-infinite 
system. A time-dependent Green's function at the boundary, for obeys the RG equation: 

[d/d]nT + /3(g)d/dg + 2 1 (g)]G(T,g) = 0, (4.1) 

where 7 is the anomalous dimension of the boundary operator whose Green's function is being calculated. This is 
given, to O(g), by Eq. (2.2). g(r) in Eq. ( [4.1| ) is the effective coupling constant at scale r. Solving this equation we 
obtain: 

(lnM)- 46 ™/'' 

< 0(r, 0)0(0, 0) >^ . (4.2) 
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Both the exponent, x n and the power of the logarithm are different than what occurs in the bulk. For the lowest 
dimension boundary operator of spin S the factor of ln|r| is raised to the power Sl(Sl + 1)- In particular, for the 
spin operator at the boundary in the lattice Heisenberg model, the correlation function behaves as: 

< S (t) • 5 (0) >^ constan t , ' U . (4.3) 

m 

The imaginary part of the retarded Green's function at zero frequency and finite T, obtained from the Fourier 



transform of Eq. (4.3), gives the nuclear magnetic resonance relaxation rate, 1/Ti for a chain with non-magnetic 



impurities. This behaves as: 

1/Ti oc T[hi(T /T)] 2 , (4.4) 

for some temperature scale To of order the exchange energy. The 1/r 2 power law was first derived in [ ^| , without 
consideration of logarithmic corrections. The discussion of the linear power- law in 1/T 1; resulting from performing 
the Fourier transform, was first discussed, as far as we know, in [ [ll|. These authors also attempted to calculate 
the logarithmic correction. However, their result, for which no derivation was given, differs from ours, containing 
[ln(T /T)] 4 rather than [lnT /T)] 2 . This behaviour is to be contrasted with that for the pure system, in which 
1/Ti cx [h^T/To)] 1 / 2 , constant up to a log correction.^ 

Of course, an actual experiment on a doped S=l/2 chain compound would presumably average over all distances 
from the chain ends. (This is related both to the fact that the relaxing nuclei can be at arbitrary locations and that 
even a nucleus near the end of a chain will have a transferred hyperfine interaction with spins further away from the 
chain end.) We note thai- at T = and ignoring log corrections, the spin self-correlation for a spin a distance x from 
the chain end is given byu: 

where v is the spin- wave velocity. At sufficiently long times this decays as 1/t 2 for all x. However, for t << x/v 
it exhibits the bulk behaviour, decaying as 1/t. Thus we expect that the zero frequency finite T Fourier transform, 
which determines 1/Ti, will be essentially constant (ignoring log corrections) down to a temperature of order v/x, 
below which it will vanish .essentially linearly in T. 

1/Ti has been measurecO for the quasi oae dimensional antiferromagnetic compound Sr2Cu03 obtaining apparent 
agreement with the field theory prediction,tl (InT) 1 / 2 . Broad shoulders observed in the NMR intensityEl were inter- 
preted as resulting from the distribution of local susceptibilities predicted by field theory methodsE] for chains with 
free ends. Possibly such data will also verify the distribution of relaxation rates resulting from impurities. 

After this work was completed we managed to obtain a copy of a three year old preprint J 2 ! which was never 
published nor available on the xxx archive, and which derived many of the results obtained here. 
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